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ABSTRACT 

We present Chandra gas temperature profiles at large radii for a sample of 13 nearby, relaxed galaxy clusters 
and groups, which includes A133, A262, A383, A478, A907, A1413, A1795, A1991, A2029, A2390, MKW4, 
RXJ1159+5531, and USGC S152. The sample covers a range of average temperatures from 1 to 10 keV. 
The clusters are selected from the archive or observed by us to have sufficient exposures and off-center area 
coverage to enable accurate background subtraction and reach the temperature accuracy of better than 20-30% 
at least to r = 0.4 - O.Srigo, and for the three best clusters, to 0.6 - O.Vrigo. For all clusters, we find cool gas in 
the cores, outside of which the temperature reaches a peak at r ~ 0.15 rigo and then declines to ~ 0.5 of its peak 
value at r ^ 0.5ri8(). When the profiles are scaled by the cluster average temperature (excluding cool cores) and 
the estimated virial radius, they show large scatter at small radii, but remarkable similarity at r > 0.1 - 0.2ri8o 
for all but one cluster (A2390). Our results are in good agreement with previous measurements from ASCA 
by Markevitch et al. and from Beppo-SAX by DeGrandi & Molendi. Four clusters have recent XMM-Newton 
temperature profiles, two of which agree with our results, and we discuss reasons for disagreement for the other 
two. The overall shape of temperature profiles at large radii is reproduced in recent cosmological simulations. 
Subject headings: clusters: general 



1. INTRODUCTION 

Radial temperature profiles of the hot intracluster medium 
(ICM) in galaxy clusters and groups is one of the prime tools 
to study the gravitational processes responsible for large-scale 
structure formation and non-gravitational energy input into 
the ICM. The temperature profile is also an important cos- 
mological measurement because in dynamically relaxed sys- 
tems, it is the basic ingredient in estimating the total clus- 
ter mass distribution assuming hydrostatic equilibrium of the 
ICM. X-ray mass measurements at radius r are only as accu- 
rate as T and dT/dr at that r (e.g., Sarazin 1988). 

Temperature measurements at large cluster radii are tech- 
nically challenging. Cluster brightness at the virial radius 
is only 10% of the Cosmic X-ray Background (CXB) in the 
soft X-ray band, and a smaller fraction of the total (detector 
+ CXB) background. The background spectrum is usually 
much harder than that of the cluster, leading to even lower 
surface brightness contrast at higher energies. The practical 
implication is that accurate spectral analysis near the virial ra- 
dius is not possible, even with very long exposures, unless the 
background can be subtracted with better than 1 % accuracy. 
Such accuracy is unachievable with past or present X-ray tele- 
scopes. The surface brightness contrast near 0.5 of the virial 
radius is higher, and the required accuracy of the background 
subtraction is ~ 3%. This is easily achievable with Chandra 
(Markevitch et al. 2003) and marginally feasible with XMM- 
Newton (Nevalainen, Markevitch, & Lumb 2005). 

Cluster observations with telescopes operating below the 
Earth's radiation belts, such as ASCA and Beppo-SAX, are 
less affected by the background. However, poor point spread 
functions (PSF) of AS' CA and Beppo-SAX, l'-2', was a major 
problem for temperature profile measurements. This problem 
is non-existent for Chandra and not a particular concern for 
XMM-Newton except in the cores of clusters with peaked sur- 
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face brightness profiles (e.g., Markevitch 2002). 

Given the technical difliculties, early measurements of clus- 
ter temperatures at large radii have been controversial. ASCA 
was the first instrument with the necessary spectro-imaging 
capability. Cluster analysis was complicated by the complex 
PSF of the ASCA telescope, which if not taken into account, 
resulted in radially increasing apparent temperatures. Using 
ASCA data, Markevitch et al. (1996; 1998, hereafter M98; 
1999) obtained temperature profiles for a sample of 32 nearby 
clusters, which showed significant declines with radius be- 
tween r = 0.1 - 0.6rvir (hereafter, rigo s rigo from Evrard, 
Metzler & Navarro 1996). In clusters without obvious merg- 
ers, the radial temperature profiles outside the cool cores were 
similar when normalized to the virial radius. Consistent re- 
sults were obtained, e.g., by Ikebe et al. (1997), Cannon et al. 
(1999) and Finoguenov et al. (2001). However, White (2000) 
found that most clusters in their large ASCA sample "are con- 
sistent with isothermality at the 3cr confidence level." Some 
of the differences between White (2000) and other works can 
be attributed to White's use of a PSF model that overestimated 
scattering at low energies. Also, White typically did not ex- 
tend measurements to large radii because of large uncertainty 
inherent in his image deconvolution method, and in the over- 
lapping radial range most of his temperature profiles are in 
fact consistent with M98. 

Beppo-SAX was another instrument capable of spatially re- 
solved spectroscopy; compared to ASCA, it had better PSF. 
Irwin & Bregman (2000) analyzed 1 1 clusters and reported 
isothermal or even increasing radial temperature profiles. 
However, De Grandi & Molendi (2002) pointed out a tech- 
nical error in that work. Instead, De Grandi & Molendi 
found, in their analysis of 21 clusters, declining tempera- 
ture profiles, in good agreement with M98 outside the cores 
(r > 0.15 - 0.2 rvh), although less peaked closer to the cen- 
ter Similar profiles using Beppo-SAX data were obtained by 
Ettori et al. (2000) and Nevalainen et al. (2001). 

Temperature profiles for individual clusters have been de- 
rived in the past few years using both Chandra and XMM- 
Newton observations. Previous Chandra results were mostly 
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TABLE 1 

Cluster sample and Chandra OBSERvmoNS 



Chandra observations 



Cluster 




z 


Nh° 


Aim point 


Exposure, ksec 


ACIS Mode 


'"max /n 80 


Comments 


A133 


4 2 


0.057 


15^x1 0^" 


S+I 


40+90 


F+VF 


0.67 


1 4' oflFset in ACTS-T 


A262 


2.1 


0.016 


8.10 X lO^o (*) 


S 


30 


VF 


0.30 




A383 


4.9 


0.188 


3.92 X lO^" 


I+I+S 


20+10+20 


VF 


0.42 




A478 


7.9 


0.088 


variable (*) 


s 


43 


F 


0.60 




A907 


5.9 


0.160 


3.87 X lO^" (*) 


1 


49+35+11 


VF 


0.58 




A1413 


7.3 


0.143 


2.19X lO^o 


I 


10+113 


VF 


0.70 




A1795 


6.1 


0.062 


1.19 X 10^0 


S+I 


7x15 


VF 


0.53 


Multiple offsets 


A1991 


2.6 


0.059 


2.45 X 10^" 


s 


36 


VF 


0.48 




A2029 


8.5 


0.078 


3.04 X 10^1 


S+I 


87+9 


F+VF 


0.60 




A2390 


8.9 


0.230 


1.07 X 10^1 (*) 


s 


95 


VF 


0.90 




MKW4 


1.6 


0.020 


variable (*) 


s 


30 


VF 


0.39 




RXJ 1159+5531 


1.9 


0.081 


1.20 X 10^" 


s 


70 


VF 


0.38 




USGCS152 


0.7 


0.015 


1.55 X 10^' 


s 


30 


VF 


0.45 





^ — Emission-weighted temperature (keV), excluding central 70 kpc (§5). 

* — Galactic absorption column density (cm^^) adopted in this paper (§3). Stars mark those clusters with different radio and X-ray values of Nh- 



confined to the central region, typically within 0.2-0.3 virial 
radii (David et al. 2001; Schmidt, Allen & Fabian 2001; Maz- 
zotta et al. 2002; Sun et al. 2003a; Johnstone et al. 2002; 
Lewis, Buote & Stocke 2003; Buote & Lewis 2004; Sun 
et al. 2003b) although attempts were made to go to larger 
radii (Markevitch et al. 2000b; Markevitch & Vikhlinin 2001). 
These measurements cannot be used to test the temperature 
decline at large radii observed by ASCA and Beppo-SAX. Pub- 
lished XMM-Newton temperature profiles extend to a larger 
fraction of the virial radius. Declining temperature profiles 
are observed in some clusters (Pratt & Arnaud 2002; Taka- 
hashi & Yamashita 2003; Zhang et al. 2004; Piffaretti et al. 
2004; Belsole et al. 2005). However, there are also clusters 
reported to be isothermal at large radii (Majerowicz, Neu- 
mann & Reiprich 2002; Pointecouteau et al. 2004; Belsole et 
al. 2004; Sakelliou & Ponman 2004; Pratt & Arnaud 2004). 
These lists include results for both relaxed clusters and major 
mergers, so some of the differences could be explained by the 
differences in the cluster dynamical state. 

Chandra is well-suited for measurements of the tempera- 
ture profiles to 0.5-0.6 of the virial radius thanks to its sta- 
ble detector background, and fine angular resolution. In this 
paper, we present measurements of the temperature profiles 
to at least ~ 0.4 of the virial radius for a sample of 13 low- 
redshift clusters observed by Chandra with sufficiently long 
exposures that the statistical temperature uncertainties at this 
radius are reasonably small. The clusters are listed in Table 1. 
All these objects have a very regular overall X-ray morphol- 
ogy and show only weak signs of dynamical activity, if any. 
The main goal these observations was mass determination at 
large radii from the hydrostatic equilibrium equation. All of 
them have sufficient off-center area coverage to enable accu- 
rate background modeling and subtraction, which is the crit- 
ical element of our analysis. Three clusters (A133, A907, 
A1413) were observed by us with a specific setup optimized 
for studying the outer regions. Even though the present cluster 
sample is not unbiased, it represents an essential step towards 
reliable measurements of the gas density and temperature pro- 
files to a large fraction of the virial radius. 

We assume h = 0.72, = 0.3, Q.\ = 0.7. Measurement 
uncertainties correspond to 68% CL. 



2. CHANDRA DATA ANALYSIS 

The nominal aim points (ACIS-I or ACIS-S) and expo- 
sure time for Chandra observations analyzed in this paper are 
listed in Table 1 . Most of the observations were telemetered 
in VFAINT mode which provides for better rejection of the 
particle-induced background. In the case of ACIS-I pointings, 
we used the data from four ACIS-I chips for the temperature 
analysis, and S2 chip was generally used for monitoring the 
background. In the ACIS-S pointings, we used S3, S2, and 
available ACIS-I chips. The SI chip was used mainly to mon- 
itor the background, but in low-redshift cool clusters (A262, 
MKW4, USGC S152), we used the SI data also for the tem- 
perature profiles. In the rest of this section, we outline our 
Chandra data reduction procedures. 

2. 1 . Calibration Corrections to Individual Photons 

We start with "level 1" photon lists and apply standard pro- 
cessing using the CIAO tool acis_process_events. This 
includes correction for Charge Transfer Inefficiency (CTI; 
Townsley et al. 2000, Grant et al. 2004), re-computation of 
event grades, and detection of so-called afterglow events. We 
then remove photons detected in bad CCD columns and pix- 
els and also those with bad ASCA grades (1,5,7). For ob- 
servations telemetered in VFAINT mode, we apply additional 
background screening by removing events with significantly 
positive pixels at the border of the 5x5 event island. 

The next step is to compute calibrated photon energies. We 
apply updated ACIS gain maps and correct for its time depen- 
dence. Both corrections are available in the CIAO 3.1. The 
new gain maps are required for analyzing SI data; improve- 
ments are very small for other chips. 

The final step is to examine background light curves dur- 
ing each observation in order to detect and remove the flaring 
episodes. The flare detection was performed following the 
reconomendations given in Markevitch et al. (2003). Clean 
exposure times are listed in Table 1. 

2.2. Spectral Response Calibration 

To take into account the spatial dependence of effective area 
and detector energy resolution, our analysis uses several im- 
portant recent cahbration advances outHned below. 
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Fio. 1. — Composite Chandra images in the 0.7-2 keV band. The images were not flat-fielded to indicate the diflerence in exposure coverage for diff'erent 
locations within the cluster Each CCD chip is 8' on a side. North is up and East is to the left. 



ACIS optical blocking filters have been contaminated in- 
flight by a substance containing C, O, and F. Photoelectric ab- 
sorption in the contaminant strongly reduces the low-energy 
effective area. Thanks to a number of calibration analyses, 
the properties of the contaminant are now known sufficiently 
accurately. The absorption spectrum of the contaminant has 
been measured using grating observations of bright AGNs 
with continuum spectra (Marshall et al. 2004). Time depen- 
dence of its thickness in ACIS-S was determined from the flux 
ratios of the calibration source lines. These corrections are 
built into CIAO since v3.0. New to our analysis is the inclu- 
sion of the spatial distribution of the contaminant. The center- 
to-edge difference of its optical depth is At ^ 0.5 around 
0.7 keV in bofli ACIS-I and ACIS-S arrays. Such variations, 



if unaccounted for, may have a major impact on the derived 
ICM temperatures. Recently, the spatial distribution of the 

contaminant has been measured accurately, so that the resid- 
ual variations of the effective area are less than 3% at ener- 
gies above 0.6 keV (Vikhlinin 2004a). This correction was 
released to the Chandra users in October, 2004. 

An additional non-uniformity of the CCD quantum ef- 
ficiency arises because for a certain fraction of photons, 
the trailing charge produced by CTI makes them appear as 
charged particle-induced events, and they are screened out. 
Earlier cahbration of this effect was based on a limited dataset 
and lacked an accurate model for its energy dependence. 
The resulting quantum efficiency residuals are ^ 5% near 
E — 6keV, marginally significant for the ICM temperature 
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A2390 



RXJ 1159+5531 




Fio. 1. — continued. 

measurements. The new calibration released with CIAO 3.1 
eliminates these residuals (Vikhlinin 2004b). 

The same CIAO release corrected an error in the CCD av- 
erage quantum efficiency (QE) for the backside illuminated 
chips (S3 and SI) both at low and high energies (Edgar et al. 
2004). This error, if neglected, results in a ^ 10% mismatch 
in the temperature values determined from the back- (BI) and 
fronlside illuminated (FI) chips. Temperature profiles span- 
ning both types of chips would be affected. 

Finally, the released calibration underestimates the effec- 
tive area of the Chandra mirror by ~ 10% just above the Ir M 
edge, probably because the mirror surface is contaminated by 
a thin hydro-carbon layer (Marshall et al. 2003). We find that 
the corresponding effective area correction can be approxi- 
mated by a "positive absorption edge", 

C = exp(0.15(£/£o)"^) for E > Eq, (1) 

where £0 = 2.07 ke V. This correction reduces the best-fit tem- 
peratures by 5%. The expected effect of mirror contamination 
on the vignetting is negUgible (D. Jerius, private communica- 
tion). Corrections of the mirror area and updates to the QE of 
BI chips supersede the fudge factor of 0.93 that was suggested 
for the FI QE (Markevitch & Vikhlinin 2001). 

With all these recent caUbration updates applied, the vari- 
ations of the effective area within the Chandra field of view 
are modeled accurately to within 3%. This is verified directly 
by comparison of the A 1795 spectra observed at different lo- 
cations (Vikhlinin 2004a), and indirectly by the agreement of 
the temperature profiles derived in different CCD chips. 

2.3. Background Subtraction 



MKW4 



USGC SI 52 




As was discussed above, correct background subtraction 
is crucial for temperature measurements at large radii. The 
ACIS detector background can be modeled very accurately, 
to within 2% rms, which is one of the major advantages of 
Chandra over XMM-Newton for studying cluster outskirts. 

The baseline background model can be obtained by using a 
compilation of the blank-field observations, processed identi- 
cally to the cluster data, and "re-projected" onto the sky using 
the aspect information from the cluster pointing (see Marke- 
vitch et al. 2000 for a detailed description). This procedure 
results in background residuals of < 10% - 20% outside the 
obvious flaring periods. The accuracy can be significantly im- 
proved by small adjustments to the baseUne model. 

2.3.1. Quiescent Background 

ACIS background above 2 keV is dominated by events from 
the charged particles. There are secular and short-term varia- 
tions of the intensity of this component by as much as 30%, 
but its spectrum is very stable. Therefore, these trends can 
be accounted for by simply changing the background normal- 
ization. The renormalization factor can be derived for each 
observation using the data in the 9.5-12 keV band where the 
Chandra effective area is nearly zero and so all the observed 
flux is due to the background. Such a renormalization reduces 
the uncertainties in the modeling of the quiescent background 
to 2% rms, provided that flaring periods were screened out 
properly (Markevitch et al. 2003). 

We extract and fit cluster spectra separately in the BI and 
FI chips, and flare filtering is also specific for each CCD type. 
Therefore, the background renormalizations were determined 
independently for these chip sets. A 2% scatter in the overall 
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normalization per set per pointing was included in the overall 
temperature uncertainty budget. The accuracy of the quies- 
cent background subtraction was verified using the spectra in 
the regions far from the cluster centers. 

2.3.2. Soft Diffuse X-ray Background 

In addition to the particle-induced background, the blank- 
field datasets contain the diffuse X-ray background. This 
component makes a large contribution to the total background 
below ~ 1 keV. The soft X-ray background was studied with 
Chandra by Markevitch et al. (2003) and XMM-Newton by 
Lumb et al. (2002). For most of the sky above the Galac- 
tic plane, the background spectrum is well represented by the 
MEKAL model with the Solar metallicity and T » 0.2 keV. 
A second component with T » 0.4 keV appears in locations 
where the diffuse flux is high. A fraction of the soft back- 
ground is geocoronal in origin and variable in time (Wargelin 
et al. 2004). This component is dominated by the O VII and 

VIII lines and its spectrum also can be approximated with 
sufficient accuracy by a thermal plasma model. 

The blank-field datasets contain a typical mixture of the 
Galactic and geocoronal backgrounds. Since these compo- 
nents vary with location or time, the soft background in in- 
dividual pointings is usually slightly different. Therefore, 
appropriate adjustments are needed, even though they rarely 
lead to qualitative changes in the derived temperature profile. 

A certain fraction of the detector area in most of our obser- 
vations was essentially free from the cluster emission and so 
we were able to determine the soft background adjustments 
in situ. We extracted spectra in those regions, subtracted the 
adjusted blank-field background (§2.3.1), and fit the residu- 
als with an unabsorbed MEKAL model, whose normalization 
was allowed to be negative. The fit was performed in the 0.4- 

1 keV band because the soft background is usually dominated 
by the oxygen lines near 0.6 keV (Fig. 2). The derived ad- 
justment is the real sky X-ray emission and is subject to the 
spatial variations in the effective area. Therefore, it must be 
included as an additional component in the spectral fits, with 
its normalization scaled by the region area. This approach 
was used before by, e.g., Markevitch & Vikhlinin (2001), and 
a similar method is used for XMM-Newton (e.g., Pratt et al. 
2001, Majerowicz et al. 2002). 

The soft-band adjustments are within 20% of the nominal 
background flux in the 0.6-1 keV band in all our clusters, and 
do not qualitatively change the behavior of the temperature 
profiles in the radial range we consider. The only exceptions 
are A1991 and A2029 which are projected on the North Polar 
Spur. The soft background, if neglected, would have a serious 
impact on the derived temperature profiles in these cases. The 
statistical uncertainty in the normalization of this component 
is included in the uncertainties for the cluster temperatures'*. 

Our adjustment procedure implicitly assumes that there are 
no significant fluctuations of the soft background on ~ 15' 
scales. We cannot test this assumption internally in each ob- 
servation because cluster emission fills large fraction of the 
field of view. This is a limitation which can be addressed only 
by accurate mapping of the soft background in the vicinity of 
each cluster. However, we do not expect this to be a serious 
problem because the adjustments are usually small. In those 

We fitted the cluster spectra with the background normahzations fixed 
at ±1(T from the nominal value. The changes in the best-fit cluster T were 
treated as the corresponding uncertainty. They were added to the statistical 
uncertainties in quadrature. Identical technique was used to estimate the un- 
certainties related to normalization of the quiescent background. 




0.5 1 

E, keV 



Fig. 2. — The spectrum of the oft'set A133 pointing in the S2 chip with 
the nominal background subtracted. Residuals give an example of oversub- 
traction of the sky soft background. The fit is MEKAL model with negative 
normalization, Solar metallicity, and T = 0.18 keV. It is included in the final 
background model. 

cases with the high Galactic foreground flux, we have used the 
ROSAT All-Sky Survey maps to verify the absence of surface 
brightness gradients on relevant spatial scales. 

2.3.3. Subtraction of Readout Artifact 

ACIS CCDs are exposed to the source during the readout 
cycle, which takes 41 ms or 1.3% of the nominal exposure 
time (3.2 s). This results in 1.3% of the source flux being 
uniformly re-distributed into a strip spanning the entire CCD 
along the readout direction, creating a "readout artifact". This 
effect must be accounted when the readout artifact from the 
bright, central region contaminates the low surface bright- 
ness cluster outskirts. Fortunately, the readout artifact for 
sources without pile-up can be subtracted almost precisely, 
using a technique proposed by Markevitch et al. (2000a). A 
new dataset is generated with the CHIPY-coordinate of all the 
photons in the original observation randomized, and their sky 
coordinates and energies recalculated as if it were a normal 
observation. The obtained dataset is normalized by the ratio 
41 ms/3.2 s and treated as an additional background observa- 
tion. The normalization of the blank-field background is re- 
duced by the same 1.3% to account for this subtraction. This 
correction was applied for all our clusters. 

2.4. Spectral Fitting 

The cluster spectra were extracted in annuli centered on the 
X-ray surface brightness peak, separately for each pointing 
and for FI and BI CCDs. The radial boundaries were chosen 
so that rout/Hn = 1.5 except for A262 where the statistical 
accuracy was sufficient to measure temperatures in narrower 
annuli, rout/Hn = 1-33. The spectra were fit in the 0.6-10 keV 
band to a single-temperature MEKAL model. The metallicity 
was allowed to be free in each annulus in the central region, 
where the statistical uncertainties were AZ < 0.2. At larger 
radii, metallicity was fixed at the average value in the last two 
radial bins where it was actually measured. If several datasets 
(pointings or FI/BI chips) contributed to the same annulus. 
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Fig. 3. — Projected temperature and metallicity profiles for A133. The vertical bars are at the emission-weighted radius within each annulus. The independent 
temperatures measurements from S3 and FI chips in the ACIS-S pointing, and from ACIS-I pointing are shown in green, magenta, and blue, respectively. The 
values from the joint fit to these three datasets are shown in black. The metallicity at r > 600 kpc was fixed at the average value from r = 270 - 600 kpc, Z = 0.2 1 , 
and hence shown without the vertical eiTor bars. 



they were fit jointly, with the spectral parameters (tempera- 
ture, metallicity. Galactic absorption) tied and normalizations 
free for each dataset. 

Absorption was fixed unless a significant variation was 
found between the best-fit values in three cluster-centric an- 
nuU with radii 0-200, 200-400, and 400-800 kpc, where it 
usually could be measured accurately. Significant variations 
were detected only in A478 and MKW4. Furthermore, in 
most cases, the best-fit value of Nh from the X-ray fit was 
consistent with the value from the radio surveys (Dickey & 
Lockman 1990), and we used the latter in those cases. 

3. RESULTS FOR INDIVIDUAL CLUSTERS 

3.L Abelll33 

Abell 133 is a relaxed cluster at z = 0.057 with the aver- 
age temperature T ^ 4.5 keV. It was observed by Chandra 
for 40 ksec in ACIS-S and for 90 ksec in the offset pointing 
in ACIS-I, which we have designed specifically for this mea- 
surement (Fig. 1). The X-ray image of A133 is nearly axially 
symmetric on large scales. The only detectable substructures 
are located within the central 40 kpc and likely caused by ac- 
tivity of the central AGN (Fujita et al. 2002) 

The quiescent background had to be adjusted by +5% in 
the ACIS-I pointings, any by -11% and -8% in the ACIS-S 
pointing in the FI chips and BI chips, respectively. The spec- 
trum extracted outside 20' of the cluster center shows signif- 
icant negative residuals around 0.6-0.7 keV, indicating that 
the soft X-ray background is over-subtracted. The residu- 
als can be fit by the MEKAL model with T = 0.18 keV 
and normalization corresponding to the 0.7-2 keV count rate 
(-1.9 + 0.3) X 10"^ cnt s"' arcmin"^. This component, with 
the normalization scaled by the area of the spectrum extrac- 
tion region, is included in the cluster spectral fits. 

Galactic absorption shows no statistically significant de- 
viations and the cluster-average value from the X-ray fit, 
Nh - (1-6 + 0.4) X 10^" cm"^, agrees well with the radio 
value 1.53 x 10^*^ cm"^, which was used thereafter. 



The temperature and metallicity profiles are shown in 
Fig. 3. There is excellent agreement between the tempera- 
tures measured independently from the S3 and FI chips in the 
ACIS-S pointing (green and magenta, respectively), and from 
the off-set ACIS-I pointing (blue). We find similar agreement 
between different pointings and chips types in all our clusters. 

The temperature values from the joint fits are shown by 
thick black error bars. There is a significant temperature de- 
cline at large radii, from T - 4.6 + 0.15 keV at r = 170 kpc to 
r = 1.8 + 0.6 keV in the outermost annulus at r ^ 1000 kpc. 
The upper limit on the temperature in this bin is 2.7 keV at 
the 90% CL. The temperature uncertainties quoted here and 
shown in Fig. 3 include uncertainties in the normalizations of 
the high-energy and soft background components, contribut- 
ing +0.35 and +0.1 keV, respectively in the outermost bin. 

The ICM metallicity is nearly Solar in the center and de- 
creases to Z ^ 0.2 at large radii. Similar metallicity profiles 
were derived from Beppo-SAX observations of several "cool- 
ing flow" clusters (De Grandi & Molendi 2001). Chandra 
data cannot constrain the abundance outside r - 600 kpc and 
so temperatures in this region are derived with the abundance 
fixed at Z = 0.21, the average value for r - 270 - 600 kpc. 
If, instead, one assumes that the observed abundance gradi- 
ent continues at r > 600 kpc, the best-fit temperatures will 
be only ~ 3% lower The temperature variation with metal- 
licity is much smaller than the statistical uncertainties for any 
Z < 0.5, and we will ignore it for this and all other clusters. 

Inside the central 70 kpc, we observe a sharp drop in the 
projected temperature, which is most likely caused by radia- 
tive cooling. Interestingly, there is also a spike in metallic- 
ity within this radius. Such behavior is observed in all our 
clusters. However, we note that a reliable metallicity anal- 
ysis in the central regions, where the temperature gradients 
are strong, requires spectral deprojection. This is beyond the 
scope of this work, since we are primarily interested in the 
temperature profiles at large radii. Our values of the metallic- 
ity in the very centers should be treated with caution. 
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Fig. 4. — Temperature and metallicity profiles for A1413. The blue error bars in the left panel show XMM-Newton results from (Pratt & Amaud 2002), 
renormalized by +10% to account for Chandra-XMM-Newton cross-calibration (see text). 



3.2. Abell 1413 

Abell 1413 (z - 0.1429) was observed in three pointings. 
One observation (OBSID 537) was almost entirely affected 
by a background flare and it was discarded from the analysis. 
The composite image from the remaining two observations 
is shown in Fig. 1. The cluster center was placed in the 13 
chip in the shorter observation (OBSID 1661, 10 ksec). In 
the longer observation (OBSID 5003, 75 ksec), the cluster 
center was deliberately placed off-axis in the 12 chip to obtain 
better coverage for the outer regions. Since the exposure times 
are very different, our results are dominated by the data from 
OBSID 5003. The X-ray image is elongated in the North- 
South direction. However, no obvious signatures of a merger 
are observed, except for a weak edge in the surface brightness 
at 400 kpc North of the cluster center. 

We do not detect any radial variations in the Galactic ab- 
sorption, and the averaged value from the X-ray spectral fit, 
Nh - (2.4 + 0.8) X 10^° cm"^, is consistent with the radio 



Newton values by +10%. The corrected XMM-Newton tem- 
perature profile is shown by blue error bars in Fig. 4. There is 
excellent agreement between the Chandra and XMM-Newton 
measurements, except in the very center where XMM-Newton 
can be affected by poorer angular resolution. 

3.3. Abell 1795 

Abell 1795 (z = 0.0622), was observed by Chandra multi- 
ple times, within the Guaranteed Time program in 1999 and 
early 2000, and subsequently in 6 pointings as a calibration 
target. We do not use the 1999 observation (OBSID 494) be- 
cause accurate gain calibration is still unavailable for that time 
period. The remaining 7 pointings, in which the cluster was 
placed at different locations within the S3 and 13 chips, have 
a combined exposure of more than 100 ksec (Fig. 1). 

There is no significant variation of Galactic absorption with 
radius, and the cluster-average value, Nh - (1.35 + 0.12) x 



value, 2.19 x lO^o cm-^. The spectrum from r > 1700 kpc ^^^^^ ^^^'^ ^q2o ^^^-2 ^e do not detect any signatures of 



shows negative residuals near 0.6 ke V, therefore the soft back 
ground is over-subtracted. The corresponding background 
correction can be modeled with the MEKAL spectrum with 
T - 0.20 keV and normalization corresponding to the 0.7- 
2 keV count rate (-2.7 + 0.2) x 10"^ cnt s"' arcmin"^. 

The projected temperature and metallicity profiles are 
shown in Fig. 4. The results are qualitatively similar to A133 
— there is a temperature drop by a factor of ~ 2 in the out- 
ermost bin. There is also a cool region within the central 
100 kpc associated with a spike in the metallicity profile. 

A 14 13 is one of the clusters with the most accurately mea- 
sured XMM-Newton temperature profiles (Pratt & Arnaud 
2002). Comparison with our results is obviously in order. The 
average XMM-Newton temperature of the cluster is ^ 10% 
lower than our measurements. A similar difference is ob- 
served in distant clusters (Kotov & Vikhlinin, in preparation), 
and therefore most likely reflects cross-calibration problems. 
This difference is unimportant for comparison of the tem- 
perature gradients and we uniformly renormalized the XMM- 



incorrect subtraction of the soft background at the largest ra- 
dius covered by Chandra pointings. Therefore, no additional 
background correction is done, but a typical uncertainty in 
normalization of such a component (see discussion on A133 
and A1413 above) is still included in the final error budget. 

The temperature and metallicity profiles are shown in 
Fig. 5. There is a significant temperature decrease at large 
radii relative to the temperature just outside the central cool 
region. The temperature drop is weaker than that in A133 
and A1413, but we note that Chandra observations of A1795 
cover a smaller fraction of the cluster virial radius. 

Multiple observations of A 1795 at different locations on the 
detector make this cluster an ideal case for testing the calibra- 
tion uncertainties. The temperatures derived independently 
from the data in the Bl and Fl chips are shown in Fig. 5 by 
green and blue error bars, respectively. There is excellent 
agreement, even though these datasets are from CCDs with 
substantially different quantum efficiencies, and the data at 
the same distance from the cluster center sample different de- 
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Fig. 5. — Temperature and metallicity profiles for A1795. Independent T measurements from the BI and FI chips are shown by green and blue, respectively. 
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Fig. 6. — Comparison of Chandra and XMM-Newton temperature profiles 
for A 1795. The XMM-Newton results from Arnaud et al. (2001) are shown 
in grey. The same observation was also analyzed by Nevalainen et al. (2004; 
blue error bars with caps). Note that XMM-Newton values were renormalized 
by +10% to account for cross-calibration (see § 3.2). 



tector regions and therefore different contamination depths. 

Temperature profile for A 1795 was the first one published 
from XMM-Newton (Arnaud et al. 2001). XMM-Newton re- 
sults are compared with our measurements in Fig. 6. No ob- 
vious trend for the temperature to decline with radius is sug- 
gested by the Arnaud et al. analysis (grey error bars), but the 
difference with our results is within their statistical uncertain- 
ties. A subsequent reanalysis of the same XMM-Newton data 
by Nevalainen et al. (2004), who performed a more rigorous 
background flare exclusion, produced a temperature profile 
fully consistent with our results (blue error bars with caps). 



3.4. Abell262 

Abefl 262 (z = 0.0162, T ^2 keV) was observed in ACIS-S 
for 30 ksec. The X-ray image (Fig. 1) on large scales in nearly 
symmetric, but there is substructure in the very central region, 
probably related to past activity of the central radio source 
(Blanton et al. 2004). We include data from the S 1 chip in this 
observation. This CCD is not used for temperature measure- 
ments in hot clusters, because its particle-induced background 
is very high above 5 keV, but Si's performance is comparable 
to that of S3 for clusters with T <2 keV, such as A262. 

We do not find any problems with the soft background sub- 
traction in the A262 observation. The cluster brightness is 
sufficiently high everywhere in the field of view and therefore 
any reasonable background variations have negligible effect 
on the derived temperatures. The best-fit Galactic absorption 
is constant with radius within the measurement uncertainties 
but its average value, Nh = (8.1 + 0.3) x 10^° cm~^ is signif- 
icantly higher than the radio value, 5.4 x 10^" cm"^. Excess 
absorption in this cluster was also found in the ROSAT data 
(David, Jones, & Forman 1996). David et al. demonstrated 
that the excess absorption can be fully explained by interven- 
ing Galactic cirrus visible in the IRAS lOOjt/ image. We use 
the best-fit X-ray value of Nh in the further analysis. 

The temperature and metallicity profiles of A262 are shown 
in Fig. 7. Note a strong central spike of metallicity, which 
coincides with the temperature decrement within the inner 
50 kpc. There is a significant temperature decrease at r > 
100 - 200 kpc. Note that in this case, we use a very different 
portion of the ACIS energy band for the temperature determi- 
nation and so systematics should be very different from those 
in the hotter clusters. It is reassuring, therefore, that we ob- 
serve a quaUtatively similar temperature structure. 

3.5. Abelimi 

Abell 1991 is a r ^ 2.8 keV cluster at z = 0.0592 and 
its virial radius fits almost entirely within the Chandra field 
of view. A1991 was observed for 40 ksec in ACIS-S (Fig. 1). 
Sharma et al. (2004) report substructure within the central 10" 
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Fig. 7. — Temperature and metallicity profiles for A262. 
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Fig. 8. — Temperature and metallicity profile for A1991. XMM-Newton results from Pratt & Amaud (2004) are shown in blue. Note that XMM-Newton values 
were renormalized by +10% to account for cross-cahbration (see § 3.2). 



but Otherwise, the X-ray image is very symmetric. 

Absorption, derived from the X-ray spectra, is constant with 
radius and consistent with the radio value, Nh - 2.45 x 
10^° cm"^. The cluster is projected on the North Galactic Spur 
and in this field there is a strong excess flux at low energies. 
Fortunately, the spectrum and normalization of this soft back- 
ground component can be determined from the same observa- 
tion because the Chandra's, field of view extends to 20' from 
the cluster center, where the object surface brightness is negli- 
gible. Galactic foreground was measured independently from 
the 0.4-1 .5 keV spectra in the 9.5'- 14.5' and 14.5' -20' radial 
ranges extracted in the SI and ACIS-I chips. The excess flux 
requires two components with different T, similar to results 



for the Chandra blank fields in the Spur area (Markevitch et al. 
2003). The best-fit temperatures are Ty = 0.157 + 0.029 keV 
and T2 - 0.51 + 0.15 keV, and relative emission measures 
are K2/K1 - 0.285. The normalization corresponds to a flux 
of 2.45 X 10"'^ erg s"' cm"^ per arcsec^ in the 0.4-1.5 keV 
energy band. The normalizations are consistent in all four 
spectra, indicating that this soft component indeed represents 
a uniform foreground and is not related to the cluster 

The temperature and metallicity profiles of A 1991 are 
shown in Fig. 8. The spectrum in the innermost bin clearly 
requires several temperature components and therefore the de- 
rived metallicity for a single-temperature fit is artificially low. 

XMM-Newton results for A 1991 were recently published by 
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Fig. 9. — Temperature and metallicity profiles for A383. 
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Fig. 10. — Temperature and metallicity profiles for A478. The XMM-Newton measurements (renormalized by +10%, see § 3.2) are shown in blue. 



Pratt & Arnaud (2004). They derive an essentially isothermal 
profile, which disagrees with our results at r > 300 kpc at 
a ~ 3cr confidence level (blue points in Fig. 8). The differ- 
ence in modeling of the excess soft background in this field 
and background flares are obvious suspects for the discrep- 
ancy between the XMM-Newton and Chandra results; unfor- 
tunately, Pratt & Arnaud do not provide information neces- 
sary for a more detailed comparison. 

3.6. Abell383 

Abell 383 {z - 0.1883) was observed with Chandra in 3 
pointings, for 20 and 10 ksec in AClS-1, and for 20 ksec in 
ACIS-S (Fig. 1). There are faint extended X-ray sources in 



the X-ray surface brightness at 5.5' to the South-East and at ~ 
7' to the North- West, probably related to infalling subgroups. 
They were excluded from the radial profile analysis. 

Galactic absorption derived from the X-ray spectra is con- 
stant with radius and consistent with the radio value, 3.92 x 
10^*^ cm"^. The spectrum extracted in the outermost regions 
of the field of view shows negative residuals near 0.6 keV, 
similar to those in the A133 observation (Fig. 2). The resid- 
uals can be fit with the MEKAL model with T - 0.18 keV 
and a negative normalization corresponding to a 0.7-2 keV 
count rate of (-1.5 ± 0.3) x 10"'' cnt s"' arcmin"^. This com- 
ponent was included in the spectral fits. The temperature and 
abundance profiles for A383 are shown in Fig. 9. 
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Fig. 11. — Temperature and metallicity profiles for A907. 

3.7. Abell478 

Abell 478, one of the highest-temperature clusters in our 
sample, was observed in ACIS-S for 43 ksec (Fig. 1). In the 
central 30" there is substructure associated with activity of the 
central AGN but otherwise, the cluster is symmetric. Chandra 
analysis of the central 4' was presented by Sun et al. (2003a). 
Our analysis supersedes the temperature profile in that paper 
because we use updated calibrations. This cluster was also 
studied by XMM-Newton (Pointecouteau et al. 2004). 

The spectra from the outermost regions of the field of view, 
which are almost free from cluster emission, do not show any 
significant residuals in the soft band. However, there are sig- 
nificant variations in the Galactic absorption. The best-fit Nh 
changes linearly with radius from (3.09 + 0.09) x 10^' cm"^ 
at r = to (2.70 + 0.06) x 10^' cm"- at r = 4'; these are sig- 
nificantly in excess of the radio value, 1.51 x 10^' cm"^. At 
r > 5', Nh drops quickly, reaching 1.5 x 10^' cm"^ at r 8' 
and stays constant at larger radii. Pointecouteau et al. obtained 
very similar results from the XMM-Newton data and argued 
that the excess absorption is most likely caused by a Galac- 
tic dust cloud projected on one side of the cluster, including 
the center. Since absorption is a strong function of radius in 
this case, we leave it free in all spectral fits. The absorption 
distribution may not be azimuthally symmetric relative to the 
cluster center and using an average value in an annulus is not 
accurate. We ignore this problem because Nh variations are 
slow inside the central 5' and at larger radii, the ACIS field of 
view covers only the South- West sector (Fig. 1). 

The temperature and metallicity profiles are shown in 
Fig. 10. The abundance is unconstrained by the Chandra data 
in the individual annuli at r > 500 kpc and we fixed it at the 
average value in the r - 500 - 1600 kpc range, Z - 0.33. Any 
reasonable variations of metallicity do not affect the temper- 
ature determination in hot clusters such as A478. As in other 
clusters, we observe a strong central temperature decrement, 
a peak at r ~ 200 - 300 kpc where T ^ 8.5 keV, and a gradual 
decUne to T ^ 6 keV at r > 700 kpc. The XMM-Newton tem- 
perature profile is shown in blue in the left panel of Fig. 10. It 
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is qualitatively different from Chandra profile — the temper- 
ature gradually increases with radius to r ~ 4' and then stays 
approximately constant. The biggest difference between the 
two instruments, however, is within the central 3'. Chandra 
results suggest a compact (r < 1') cold central region sur- 
rounded by a hotter annulus, while XMM-Newton observes 
a more distributed central temperature decrement. Outside 
the central 3', Chandra and XMM-Newton temperatures are 
in reasonable agreement. We note that A478 has one of the 
most centrally-peaked X-ray brightness profiles in our sam- 
ple. Therefore, is the prime suspect for the temperature dis- 
crepancy in the central region is a combination of incomplete 
correction for the PSF effects in the XMM-Newton analysis 
and the complex temperature structure of the central 100 kpc 
region in this cluster (Sanderson, Finoguenov & Mohr 2004). 

3.8. Abell 907 

Abell 907 (z = 0.1603) was observed in 3 separate ACIS-I 
pointings for 49, 35, and 1 1 ksec, designed specifically for 
this measurement (Fig. 1). There are no significant radial 
variations of the Galactic absorption, but the cluster-averaged 
value, Nh - (3.87 + 0.2) x 10^*^ cm"^, is significantly lower 
than that from the radio data, 5.4 x 10^*^ cm"^. We use the 
X-ray value in our analysis. There are also negative residu- 
als near 0.6 keV in the spectra from the outermost regions. 
The corresponding background correction can be modeled 
as a MEKAL spectrum with T - 0.278 keV and normal- 
ization corresponding to the 0.7-2 keV flux (-4.8 + 1.3) x 
10"* cnt s"' arcmin"^. The derived temperature and metallic- 
ity profiles are shown in Fig. 11. 

3.9. Abell 2029 

Abell 2029 (r = 8.5 keV, z = 0.078) was observed twice in 
ACIS-S for 10 and 80 ksec (both times in Faint mode) and in 
ACIS-I for 10 ksec (VFaint mode). We use all three pointings. 
The X-ray morphology of this cluster is very regular outside 
the central ~ 50 kpc. Chandra data for the central region of 
A2029 was used by Clarke et al. (2004) to study interaction 
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Fig. 12. — Observed spectra in the outermost region of A2029 (r > 14.5'). 
Data from S 1 and FI cliips are sliown in green and blue, respectively. The 
model includes the cluster emission and a 2-temperature Galactic component 
(see text for details). Total fit is shown by solid histograms. Dashed lines 
show contribution of the cluster component. 

between X-ray and radio-emiting plasma, and by Lewis et al. 
(2003) to derive the inner dark matter profile. 

The best-fit Galactic absorption everywhere in the central 
region is consistent with the radio value, 3.0 x 10^*^ cm"^. 
The major complication for the Chandra analysis of A2029 
is high Galactic soft X-ray foreground (A2029 is projected 
on the Northern Galactic Spur). Because of the cluster prox- 
imity, its emission is non-negligible everywhere in the ACIS 
field of view. Therefore, we had to fit the Galactic flux jointly 
with the cluster emission using the spectrum integrated at 
r > 14.5' from the cluster center. The cluster and Galac- 
tic flux are easy to separate spectrally because the latter is 
much softer (Fig. 12). We find that a 2-temperature MEKAL 
model is required to represent the Galactic foreground. The 
situation is similar to that in A 1991, another cluster in our 
sample projected on the Northern Galactic Spur Its spectrum 
and normalization is determined from the joint fit to spectra 
extracted in the SI and Fl chips. The model includes an ab- 
sorbed MEKAL spectrum to represent the cluster, and two 
additional unabsorbed MEKAL components with the Solar 
metallicity to represent the foreground. Spectral parameters 
for all components are tied for SI and Fl datasets; the nor- 
malizations of the cluster component are independent, while 
those of the Galactic components are tied to be proportional to 
the area covered. The best fit provides an excellent description 
of both datasets (solid lines in Fig. 12). The best-fit temper- 
atures of the Galactic components are Ti = 0.22 ± 0.05 keV, 
and T2 = 0.46 + 0.07 keV, and the ratio of emission mea- 
sures is K2/KI - 0.53. These spectral parameters are similar 
to those derived for A 1991 and also for Chandra blank field 
in the Spur region (Markevitch et al. 2003), which gives us 
confidence that the Galactic flux has been measured correctly. 
The best-fit temperature of the cluster component in this re- 
gion is r = 3.4;^') g keV (including systematic uncertainties of 
the quiescent detector background). 

Temperature fits in other regions are performed with nor- 
malization of the Galactic component scaled by area. The 



temperature and abundance profiles are shown in Fig. 13. 
Despite difficulties with the Galactic foreground modeling, 
A2029 has the most accurately measured temperature profile 
among hot, T > 8 keV, clusters in our sample. 

3.10. Abell2390 

Abell 2390, the hottest and highest-redshift cluster in our 
sample iz - 0.2304) was observed in 3 ACIS-S pointings 
for 10, 10, and 100 ksec. We discard the first two short 
observations. One of them was performed in 1999 and ac- 
curate gain calibration for this period is unavailable. The 
second short observation was telemetered in FAINT mode, 
which results is poorer background rejection. The X-ray im- 
age shows some substructure in the central 1 .6' or ~ 300 kpc 
(cavities and elongations) but is more regular on larger scales 
(Fig. 1). The best-fit absorption is constant with radius, Nh - 
(1.07 + 0.05) X 10^' cm"^, but is higher than the radio value 
6.8 X 10^" cm"^; we use the X-ray value in the further anal- 
ysis. The spectra from the outermost regions show positive 
soft-band residuals which can be fit with the MEKAL spec- 
trum with T - 0.24 keV and normalization corresponding to 
the 0.7-2 keV flux (8.9+2.7)x lO"*" cnt s"' arcmin"^; this was 
used as an additional soft background correction. The temper- 
ature and metallicity profiles are shown in Fig. 14. There is a 
marginally significant temperature drop from T ^ 9.5 keV at 
r ~ 500 kpc to 4.7 + 2.3 keV at r > 1200 kpc. 

3.11. RXJ 1159+5531 

RXJ 1159+5531 (T = 1.9 keV, z = 0.081) is an "X-ray 
Over-Luminous Elliptical Galaxy" selected from the 160 deg^ 
ROSAT PSPC Survey (Vikhlinin et al. 1999). Opticafly, this 
object appears as a nearly isolated elliptical galaxy but its X- 
ray luminosity and extent is typical of poor clusters. It was 
observed in AClS-1 for 20 ksec in the Fall of 1999 and in 
ACIS-S for 80 ksec in 2004. Accurate calibration is still un- 
available for the 1999 data and so we do not use this pointing 
in the present analysis. The optical structure suggests that 
this object is probably dynamically old and fully relaxed. The 
only detectable substructure in the X-ray image is, indeed, 
confined to the very central region of the cluster, r < 5 kpc. 
The average X-ray absorption is consistent with the radio col- 
umn density, Nh = 1.2 x 10^" cm"^. The soft background is 
clearly over-subtracted. The required correction is well-fit by 
a MEKAL spectrum with T = 0.17 keV and negative normal- 
ization similar to most of other similar cases. The tempera- 
ture and metallicity profiles of RXJ 1 159+5531 are shown in 
Fig. 15. 

3.12. usee SI 52 

The next two clusters are low-redshift, low-temperature 
systems. USGC S152 (z - 0.0153) was observed in a sin- 
gle ACIS-S pointing for 30 ksec (Fig. 1). This is the lowest- 
temperature {T ~ 0.8 keV) cluster in our sample. The out- 
ermost region of the field of view (r > 16.5') is virtually 
free from the cluster emission and the spectrum extracted 
at these radii shows no problems with the background sub- 
traction in the soft band. The best-fit absorption is Nh - 
(1.55 + 0.10) X 10^^ cm"^, which is significantly higher than 
the radio value, 4.34 x 10^'^ cm"^; there are no statistically 
significant variations of absorption with radius. The tempera- 
ture and metallicity profiles are shown in Fig. 16. This is the 
only cluster in our sample without a strong central tempera- 
ture decrement, probably because the energy output from a 
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Fig. 13. — Temperature and metallicity profiles for A2029. 



1500 




1000 



1500 



kpc 




3000 




2000 



Fig. 14. — Temperature and metallicity profiles for A2390. 



recent AGN outburst or from the stellar winds of the central 
cD galaxy is non-negligible in this low-temperature system. 

3.13. MKW4 

MKW4 (z = 0.0199) was observed in a single ACIS-S 
pointing for 30 ksec (Fig. 1). Faint cluster emission is de- 
tectable to the very edge of the ACIS field of view. We find 
no indication of problems with the soft background subtrac- 
tion. However, the best-fit absorption is significantly higher 
than the radio value (1.89 x 10^" cm"^) and varies with ra- 
dius from (6.4 + 0.8) x 10^" cm^^ 



near the cluster center to 
(2.6 +0.8) X 10^" cm""^ in the outermost regions. We, therefore, 
allowed Nh to be a free parameter in the spectral fits. The 



results for this cluster are shown in Fig. 17. Note a remark- 
ably strong metallicity gradient — Z changes (for a single- 
temperature model) from ~ 2 Solar in the center to < 0.3 
Solar outside the central 200 kpc. 

4. DISCUSSION OF SYSTEMATIC EFFECTS 

To the best of our knowledge, all calibration eff'ects which 
could affect the temperature profile measurements have been 
accounted for in our analysis. It is useful, nevertheless, to dis- 
cuss independent indicators of systematic errors. The main 
possible sources of bias are the following, (i) Statistical. 
Our spectral fits use minimization. The cluster emission 
is fainter at large radii and the photon counting statistics are 
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Fig. 15. — Temperature and metallicity profiles forRXJ 1159+5531 
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Fig. 16. — Temperature and metallicity profiles for USGC S152 



increasingly non-Gaussian, which in principle can lead to un- 
derestimation of the temperature. The statistical bias can be 
avoided by adequately grouping the spectra, so that there is a 
sufficient number of counts in each bin prior to background 
subtraction. The binning we used — into 26 channels of ap- 
proximately equal log width — resulted in > 100 counts per 
channel in all cases. Therefore, statistics are fully adequate 
for spectral fitting; this was verified in several cases by fitting 
the spectra with a factor of 2 coarser binning. 

(ii) Effective area. Miscalibration of the off-axis effective 
area can result in artificial temperature gradients. However, 
effective area calibration can be verified internally from our 
data. The most sensitive test is provided by different point- 



ings of A1795, with the cluster placed at different detector 
locations and off-axis angles in both ACIS-S and ACIS-I. If 
miscalibration of effective area were responsible for an ob- 
served factor of ~ 2 temperature decrease at large radii we 
would derive very different temperatures for the central re- 
gion of A1795 from on-axis and off-axis pointings. Instead, 
we find excellent agreement between these pointings (see §3.3 
and full report in Vikhlinin et al. 2004a^). Further evidence 
against effective area-related bias is provided by the off-axis 
AClS-1 observation of A133 which resulted in one of our 
most accurate temperature profiles (Fig. 3). Unlike most other 

^ ACIS calibration memo, http://cxc.harvard.edu/contrib/alexey/contmap.pdf 
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cases, the cluster center in this pointing was placed off-axis 
and the temperature decrease is observed on-axis. Finally, 
miscalibration of effective area would tend to cause all tem- 
perature profiles to be similar when plotted as a function of 
angular distance from the center Instead, we observe a very 
large scatter in such a plot (Fig. 18a). 

(iii) Background. Incorrect background subtraction can af- 
fect temperature measurements at large radii where the sur- 
face brightness is low. If incorrect background subtraction is 
the cause for temperature decrements at large radii, the back- 
ground normalization must be overestimated in all cases we 
analyzed. This is unlikely given that we used a procedure 
that was demonstrated to work very well for a large number 
of Chandra observations (Markevitch et al. 2003). Typical 
scatter of the background normalization is included in our 
temperature uncertainties. The accuracy of the background 
subtraction can be verified directly for several clusters with 
small angular size that we observed in ACIS-I, such as A383 
and A907, because a sufficient fraction of the ACIS field of 
view is free from cluster emission in these cases. There are 
no background residuals in these areas, and the temperature 
profiles for these two clusters are typical of the entire sample. 

To summarize, we cannot identify any obvious source of 
systematic error which could be responsible for the observed 
temperature decrease at large radii. An implicit indication that 
our measurements are correct is given by the similarity of all 
temperature profiles when expressed in units of the virial ra- 
dius, as discussed below. 

5. SELF-SIMILARITY OF TEMPERATURE PROFILES 

Theory predicts that clusters should be approximately self- 
similar because they form from scale-free density perturba- 
tions, and their dynamics are governed by the scale-free grav- 
itational force. Self-similarity implies, in particular, that clus- 
ter temperature (and density etc.) profiles should be similar 
when radii are scaled to the cluster virial radius, which can 
be estimated from the average temperature, rigo oc (T)'^^ (for 
detailed discussion see, e.g., Bryan & Norman 1998). This 
prediction is strongly confirmed by our measurements. 



The scaled temperature profiles for all clusters are shown 
in Fig. 18b. The temperatures were scaled to the integrated 
emission-weighted temperature, excluding the central 70 kpc 
region usually affected by radiative cooling. The same aver- 
age temperature was used to estimate the cluster virial radius, 
ri8o using a relation from Evrard, Metzler & Navarro (1996), 
ri8o = 1.95/i-'Mpc((r)/10keV)i/2. The scaled profiles are 
almost identical at r > 0.15rigo and the general trend can be 
represented with the following functional form. 



T/{T} = 



1.07, 0.035 < r/ri8o < 0.125 

1.22 - 1.2r/ri8o, 0.125 < r/ri8o < 0.6, 



(2) 



with a 15-20% scatter The only significant outlier is A2390 
(shown in magenta). However, this cluster is unusual in that 
its central cool region in this cluster extends to r ~ 400 kpc, 
probably because the cold gas is pushed out from the center 
by radio lobes. 

The strongest scatter of the temperature profiles is observed 
in the central cooling regions. This is not unexpected, be- 
cause in these regions, non-gravitational processes such as ra- 
diative cooling and energy output from the central AGNs are 
important, thus breaking self-similarity. The largest outliers 
in the centi-al region are MKW 4 and RXJ 1 159+5531 whose 
cooling regions are very compact and the temperature pro- 
files peak near r ^ 50 kpc (thus our fixed exclusion radius of 
70 kpc is too big). 

Some of the previous studies of large cluster samples with 
ASCA and Beppo-SAX have already uncovered the similar- 
ity of the cluster temperature profiles at large radii (Marke- 
vitch et al. 1998; De Grandi & Molendi 2002). Our mea- 
surements are consistent with these earlier results both qual- 
itatively and quantitatively. Red band in Fig. 19 represents 
the typical scatter of individual profiles in the M98 sample. 
Our radial trend at r > 0.1 riso is in good agreement with 
their profile, even though our sample has only two clusters in 
common with M98 (A478 and A1795). A difference in the 
qualitative trend around 0.1 ri8o is fully expected, considering 
that the ASCA measurement was at the limit of that instru- 
ment's angular resolution, and given the difference in spectral 
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Fig. 18. — Temperature profiles for all clusters plotted as a function of angular distance from the center and in units of the cluster virial radius. The temperatures 
were scaled to the cluster emission-weighted temperature excluding the central 70 kpc regions affected by cooling. The virial radii were estimated from these 
average temperatures using a relation from Evrard et al. (1996), rigg = 2.74 Mpc ((r>/10 keV)''^. Our best measurements, for A1991 (T = 2.6 keV), A133 
(T = 4.2 keV), A1413 (T = 7.3 keV), and A2029 (T = 8.5 keV), are open circles, filled circles, open squares, and filled squares, respectively. The strongest 
outliers in the right panel are MKW 4 (brown) and RXJ 1 159+553 1 (red) whose temperature profile peak at r < 70 kpc, and A2390 (magenta) whose central cool 
region extends to r ~ 400 kpc. 



modeling (we use single-temperature fits whereas M98 used 
a then-common model with ambient -t- cooling flow compo- 
nents for the cluster centers, and included only the ambient 
hot phase in the final profiles). 

The average Beppo-SAX profile for the subsample of clus- 
ters with cool cores (such as ours) from De Grandi & Molendi 
(2002) is shown by blue line in Fig. 19. The results are re- 
markably close at r > 0.2 r^Q. At smaller radii, the Beppo- 
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Fig. 19. — Compaiison of the Chandra temperature profiles and the average 
profiles from A5CA (Markevitch et al. 1998) and Beppo-SAX (De Grandi & 
Molendi 2002). ASCA results ai'e shown as the red band with the width equal 
to the scatter of the best-fit values in a sample of 30 clusters. Beppo-SAX 
results (blue line) represent the average temperature profile in a sample of 11 
"coohng flow" clusters. 



SAX profile does not recover the continuation of a tempera- 
ture increase towards smaller radii, which is consistent with 
the poorer angular resolution of that telescope. 

A declining temperature profile has been recently reported 
from an XMM analysis of 16 low-redshift clusters (Piffaretti 
et al. 2004). On average, temperature drops by ~ 30% from its 
peak value near r - 0.4ri8o (the outer radius for most clusters 
in Piffaretti et al. sample), fully consistent with our results. 

Declining temperature profiles are generally reproduced 
in recent cosmological numerical simulations (Loken et al. 
2002; Ascasibar et al. 2003; Valdarnini 2003; Borgani et al. 
2004; Ettori et al. 2004; Kay et al. 2004; Motl et al. 2004), 
although earlier simulations were more discrepant and pro- 
duced both nearly isothermal profiles (Evrard et al. 1996) 
and strong temperature gradients (Katz & White 1993). The 
best agreement with observed temperature profile seems to be 
reached by Eulerian codes (Loken et al. 2002 and discussion 
therein) and by the entropy-conserving versions of the SPH 
code (Springel & Hernquist 2002) used, e.g., by Ascasibar et 
al. and Borgani et al. A more detailed comparison with ob- 
servations requires folding the simulation output through the 
detector response (Mazzotta et al. 2004). This is beyond the 
scope of this paper. The general agreement of numerical and 
observational results suggests that the declining temperature 
profiles is a natural product of gravitational heating of ICM 
in the process of cluster formation. This conclusion is un- 
affected by the presence of heat conductivity of up to 1/3 of 
Spitzer value (Dolag et al. 2004). 

6. SUMMARY AND CONCLUSIONS 

Chandra is well-suited for measurements of the tempera- 
ture profiles to 0.5-0.6 of the virial radius thanks to its stable 
detector background, and fine angular resolution. Our analy- 
sis of high-quality observations of 1 1 low-redshift clusters re- 
veals an almost universal behavior of temperature profiles in 
the systems spanning the range of average temperatures from 
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1 to 10 keV. Projected temperature as a function of radius 
reaches a maximum at r ~ 0.1-0.2ri8(),just outside the central 
cooling region. The temperature gradually declines at larger 
radii, reaching 0.5 of its peak value near r ~ 0.5 - 0.6ri8o, the 
outer boundary of our measurements in most cases. 

In the future papers, we will use these Chandra observa- 
tions for accurate measurements of the baryon and total mass 
profiles at large radii. We will also carry out a consistent 



comparison of clusters from high-resolution numerical sim- 
ulations with their real-world counterparts. 



We thank A. Kravtsov and P. Mazzotta for useful discus- 
sion. This work was supported by NASA grant NAG5-9217 
and contract NAS8-39073. 
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